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The  strength  of  chemical  bonds  in  solids  and  liquids 
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1  Air  Force  Research  Laboratory,  Materials  and  Manufacturing  Directorate,  Wright- 
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Present  address:  Harvard-MIT  Division  of  Health  Sciences,  Massachusetts  Institute  of 
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The  strengths  of  chemical  bonds  between  atoms  are  accurately  measured  and 
widely  available  for  molecular  gases,  but  it  is  remarkable  that  a  method  to 
quantify  bond  strengths  in  liquids  and  solids  is  not  available  and  that  the  strengths 
of  these  bonds  are  generally  unknown.  We  propose  a  new  term,  the  condensed 
bond  enthalpy  (CBE),  to  specify  the  energy  contained  in  bonds  between  atoms  in 
condensed  states.  We  develop  an  approach  to  quantify  these  bond  strengths  using 
bulk  thermodynamic  and  crystallographic  data,  and  apply  it  to  generate  a  nearly 
complete  set  of  elemental  CBEs  and  a  selection  of  CBEs  between  unlike  metal  atom 
pairs.  We  demonstrate  the  validity  and  utility  of  these  values  by  applying  them  to 
several  physical  problems.  The  values  reported  here  show  a  good  predictive 
capability,  and  CBEs  from  this  approach  may  give  new  insights  into  solution 
thermodynamics  and  emerging  problems  in  the  physical  sciences. 
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Current  Bond  Strength  Representations 

The  strengths  of  chemical  bonds  between  atoms  exerts  decisive  control  in  a  wide  range 
of  physical  phenomena,  including  the  stability  and  structure  of  matter;  the  kinetics  of 
chemical  reactions  and  the  energy  required  or  produced  in  the  fonnation  of  alloys  and 
compounds.  Bond  strengths  are  commonly  quantified  by  tenns  such  as  bond 
dissociation  energy  (BDE)  and  bond  enthalpy.  Approaches  to  accurately  measure  these 
values,  typically  to  four  significant  figures,  are  well-established  ',  and  extensive 
tabulations  are  available  for  a  wide  range  of  elements  and  compounds  2’3,4,5.  While  these 
values  have  a  strong  and  demonstrated  utility,  they  are  defined  and  measured  only  for 
species  in  the  gaseous  state,  and  are  likely  to  differ  significantly  from  single  bond 
strengths  between  atoms  in  liquids  and  solids.  A  companion  approach  to  accurately 
quantify  interatomic  bond  strengths  in  the  condensed  state  does  not  currently  exist, 
leading  to  the  condition  that  surprisingly  little  quantitative  data  is  available  for  the 
strengths  of  chemical  bonds  in  liquids  and  solids. 

The  minimum  in  the  atomic  potential  energy  vs  separation  curve  gives  a  simple 
representation  of  bond  strength.  Atomic  potentials  are  often  used  in  semi-empirical 
computational  approaches.  They  are  produced  by  fitting  to  measured  values  that  include 
elastic  properties,  defect  energies  and  structural  data,  and  so  give  an  indirect  approach 
that  can  require  extensive  auxiliary  infonnation.  Pair  potentials  with  as  few  as  two 
fitting  parameters  are  used  for  simplicity,  while  more  realistic  potentials  use  as  many  as 
14  fitting  parameters  6.  Consequently,  realistic,  high  quality  atomic  potentials  are 
difficult  to  produce  and  validate,  even  in  binary  systems.  Quantitative  comparison  of 
energies  across  systems  requires  a  unifonn  approach  for  producing  atomic  potentials, 
but  fitting  techniques  are  often  unique  so  that  caution  must  be  used  in  such 
comparisons.  Ab  initio  calculations  can  give  accurate  total  energies  in  condensed 
structures,  but  it  is  difficult  to  assign  individual  bond  strengths  to  specific  atomic  pairs. 
Both  first  principles  and  semi-empirical  calculations  can  be  rather  involved,  and  it  has 
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been  difficult  to  perform  calculations  on  enough  systems  to  establish  trends  and 
patterns,  which  are  only  now  just  beginning  to  emerge  1 .  Extensive  efforts  to  predict 
phase  equilibria  in  condensed  systems  give  phenomenological  representations  of 
thennodynamic  functions,  including  composition-dependent  enthalpies  and  atomic 
interaction  parameters  .  However,  these  approaches  have  not  been  used  to  define  and 
quantify  chemical  bond  strengths.  The  regular  solution  and  quasi-chemical  models  have 
their  foundations  in  bond  energies  between  like  and  unlike  atom  pairs  in  condensed 
substances,  Sy,  and  form  a  framework  for  the  application  and  interpretation  of  these 
values  9'10.  To  our  knowledge,  £y  values  have  not  been  quantified  using  a 
thennodynamics  approach. 

Elemental  Condensed  Bond  Enthalpies 

We  develop  here  an  approach  for  estimating  bond  enthalpies  in  condensed 
elements.  The  conversion  of  one  mole  of  pure,  condensed  element  A  to  a  non¬ 
interacting,  monatomic  gas  requires  a  change  in  enthalpy  that  is  given  as  the  difference 
between  the  enthalpies  of  the  products  and  reactants 

Af H (A,  gas )  =  I /(A,  gas )  -  II(A,  cond)  ( 1 ) 

This  change  can  be  measured  as  the  heat  of  fonnation  of  a  monatomic  gas,  which  is 
numerically  equal  to  the  heat  of  sublimation  for  solid  elements  that  convert  to  a 
monatomic  gas.  Experimental  values  are  typically  given  at  298  K  and  atmospheric 
pressure  (100  kPa)  4,5,1 1,12 .  Only  changes  in  thennodynamic  values  can  be  measured,  so 
an  enthalpy  of  H(A,gas )  =  0  is  assigned  for  convenience  to  the  non-interacting, 
monatomic  gas  as  the  reference  state  against  which  changes  are  quantified.  This 
assignment  is  consistent  with  the  reference  state  typically  used  for  atomic  potential 
functions,  and  has  the  advantage  of  enabling  distinction  between  different  elements  at 
298  K.  This  differs  from  the  ‘elemental’  standard  state  often  used  in  classical 
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thermodynamics,  H(A,cond)  =  0,  which  eliminates  the  ability  to  distinguish  between 
different  elements  in  the  condensed  state  since  all  are  assigned  a  zero  enthalpy.  We  use 
the  foundational  relation  of  the  regular  solution  and  quasi-chemical  models  9,10 

H  ~  T.i,j  £ij  Pij  (2) 

where  s,j  is  the  enthalpy  of  an  i-j  atomic  bond  and  P,j  is  the  number  of  first-neighbour  i-j 
bonds.  Applying  the  previous  equations  to  one  mole  of  pure  element  A  gives 


H(A,cond )  =  - AjH(A,gas )  =  £aa(Paa)  =  sAA(NAvZAA)/2  (3) 

where  NAv  is  Avogadro’s  constant.  The  number  of  bonds  per  mole  of  atoms,  PAA,  is 
proportional  to  half  the  coordination  number,  ZAA,  to  avoid  counting  a  bond  twice. 
Rearranging  terms  gives 


Saa  =  -2  A/f/ (A,gas)/(  NA  VZAA)  (4) 

We  use  AfH(gas )  values  assessed  from  data  in  4,5,1  1,12  (see  Supplementary 
Information,  Table  SI).  A fH(gas)  is  a  positive  energy,  so  that  sAA  is  negative.  For 
convenience,  we  report  sAA  in  units  of  eV.  We  determine  sAA  for  all  elements,  excluding 
noble  gases,  Pm,  At,  Fr  and  elements  with  atomic  number  greater  than  94  where 
required  input  data  are  unavailable  (Table  1).  Bond  enthalpies  typically  range  from 
-0.1  eV  to  -1.5  eV,  with  values  as  large  as  -4.95  eV  for  elements  where  covalent 

1  T 

bonding  dominates.  These  values  are  broadly  consistent  with  expected  bond  energies 
and  show  a  clear  dependence  on  melting  temperature,  Tm  (Figure  1).  Semi-metal  and 
non-metal  elements  are  dominated  by  covalent  bonding  and  have  coordination  numbers 
much  less  than  12,  giving  strongly  negative  sAA  values  (Figure  la).  However,  all 
elements  follow  the  same  trend  when  the  same  coordination  number-  a  space-filling 
coordination  of  12-  is  used  to  calculate  sAA  (Figure  lb).  Consistent  with  Equation  4,  the 
correlations  in  Figure  1  suggest  that  an  atom  has  a  specific  capacity  for  bonding,  given 
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by  A fH{gas),  and  that  the  strength  of  the  bond  depends  on  the  number  of  bonds  formed. 
This  observation  gives  an  approach  to  quantify  the  influence  of  atomic  structure,  by 
adjusting  bond  strengths  by  the  number  of  bonds  per  atom  in  the  structure. 

Quantitative  validation  of  these  bond  enthalpies  is  given  by  comparing  measured 
vacancy  formation  energies,  enthalpies  of  fusion  and  surface  energies  with  values 
predicted  from  sAA.  The  energy  to  form  a  monovacancy  in  a  pure  element,  A///,,,  is 
modelled  as  the  product  of  eAA  and  the  net  change  in  bonds  needed  to  form  this  defect. 
Removing  an  atom  from  a  solid  interior  breaks  six  bonds  in  face-centered  cubic  (fee) 
and  hexagonal  close-packed  (hep)  structures,  and  breaks  seven  bonds  in  body-centred 
cubic  (bee)  structures.  This  atom  is  placed  on  a  free  surface  to  maintain  mass  balance, 
and  the  number  of  new  bonds  fonned  depends  on  the  crystallographic  plane  and  the  site 
occupied.  Surface  sites  include  atomically  smooth  surfaces,  atomic  ledges,  kinks  and 
surface  vacancies.  Filling  a  surface  vacancy  forms  the  largest  number  of  new  bonds,  so 
this  site  has  the  lowest  energy  and  is  strongly  preferred.  We  count  the  number  of  new 
bonds  for  fdled  surface  vacancies  on  (100),  (110)  and  (111)  planes  in  fee  structures,  for 
(0001),  (1010),  (1120)  and  (1121)  planes  in  hep  structures,  and  on  (TOO),  (110)  and 
(112)  planes  in  bee  structures.  The  average  number  of  new  bonds  fonned  is  4  for  fee 
and  hep  structures  and  4.5  for  bee  structures.  We  apply  a  small  correction  to  account  for 
the  redistribution  of  electrons  associated  with  broken  bonds  at  defect  sites  (see 
Methods).  The  results  in  Figure  2  show  that  the  predicted  energies  are  typically  within 
experimental  error. 

As  a  second  application  of  e 'aa,  we  consider  the  enthalpy  of  fusion,  A Hm,  which 
accompanies  the  transformation  from  the  crystal  to  liquid  state  at  Tm.  This  enthalpy 
change  can  be  estimated  as 


(5) 
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where  sA a  is  the  enthalpy  of  an  A-A  bond  at  Tm  and  AZ/2  is  the  change  in  the  number 
of  bonds  per  atom  associated  with  melting.  To  obtain  eAA  ,  we  adjust  AfH(A,gas )  by  the 
enthalpy  required  to  heat  the  element  to  Tm  (see  Methods).  This  reduces  saa  by  about 
10%.  We  perform  estimates  for  fee  and  hep  metals  that  do  not  undergo  allotropic 
transformations,  using  AZ=0.7  as  measured  for  liquid  Cu  14.  A  small  adjustment  to  sAa 
due  to  charge  redistribution  at  broken  bonds  (see  Methods)  is  applied.  Predicted  values 
of  A Hm  plotted  against  experimental  values  in  Figure  3  show  good  agreement,  again 
typically  within  experimental  uncertainty. 

The  final  application  of  sAa  values  given  here  is  for  solid  surface  energies,  y. 
Surface  energies  are  represented  by  the  number  of  bonds  broken  per  unit  area  (see 
Supplementary  Infonnation)  and  the  enthalpy  of  those  broken  bonds  (see  Methods). 
Surface  energies  are  predicted  for  (100),  (1 10)  and  (1 1 1)  planes  in  fee  crystals  and  for 
(100),  (1 10)  and  (112)  planes  in  bee  elements.  For  each  element,  the  mean  of  the 
predicted  values  are  compared  with  y  for  cubic  metals  measured  by  the  zero  creep 
method  at  elevated  temperatures  15.  Predictions  for  eight  elements  show  very  good 
agreement  with  values  measured  at  1723  K  or  below  and  at  0.77  Tm  or  above  (Figure  4). 
Estimates  significantly  exceed  measured  values  for  Nb,  Mo  and  W,  where  much  higher 
measurement  temperatures,  from  2273  to  2623  K,  are  used.  Enhanced  adsorption  of 
impurity  atoms  at  these  higher  temperatures  may  contribute  to  the  poor  agreement  for 
these  measurements. 

Condensed  Bond  Enthalpies  Between  Unlike  Atom  Pairs 

Consider  the  reaction  of  x  moles  of  pure,  condensed  element  A  with  y  moles  of 
element  B  to  form  one  mole  of  the  compound  AxBy  at  standard  temperature  and  pressure. 
The  enthalpy  of  formation  can  be  measured  experimentally,  and  is 

Af H(AXBV)  =  II(AXBV)  -  xH(A, cond)  -  yII(B, cond)  (6) 
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We  expand  H(AxBy )  using  Equation  2 


H(2lxSy)  =  eAAPAAfy  +  sABP;Bx°y  +  eBBPl 


AxBy  ,  _  nAxBy 


BB 


(V) 


A  B 

where  if  y  is  the  number  of  i-j  bonds  per  mole  of  AxBy.  We  replace  H{j,cond)  with 
-A fH(i,gas)  and  rearrange  terms  to  give  the  basic  relation 


£ab  — 


-Ay B 


xDy 


{[A fH(AxBy)  -  xAfH(A,gas)  -  yAfH{B,gasj\  -  sAAP/ 


AXBy 

AA 


'  AB 


£bb  Pf 


AxBy 

BB 


} 


(8) 


This  is  modified  (see  Methods)  to  give  our  final  equation 


£ab  — 


AXB 

£AA 


y 


[ &fH(AxBy)  -  xAfH(A,gas)  -  yAfH(B,gas)]  - 

(ZBy)-4BBy(pABBBy)} 


(9) 


The  terms  sa  are  replaced  with  sax  y  to  indicate  a  small  dependence  of  like-atom  bond 

A  B  A  B 

enthalpy  on  structure,  ff  '  y  is  replaced  with  y ,  the  number  of  i-j  bonds  per  unit 

cell  of  AxBy.  The  approach  for  counting  bonds  per  unit  cell  is  given  in  Methods.  U  is  the 
number  of  atoms  per  unit  cell  in  the  AxBy  structure.  The  experimentally  measured 
quantities,  AjH(AxBy)  l6’l7‘18'19  and  AjH(i,gas )  4,5‘ 1 '  ' 1 2,  are  converted  to  units  of  eV  for 
convenience. 


We  determine  Sab  for  eleven  binary  systems  where  available  data  give  a  minimum 
of  four  values  for  each  system  to  establish  compositional  trends.  We  include  data  for 
four  additional  binary  systems  with  fewer  than  four  values,  which  are  used  to  predict 
heats  of  formation  for  ternary  intermetallic  compounds  from  sAB.  We  focus  on  binary 
metallic  systems  in  response  to  Pauling’s  complaint  that  “the  great  field  of  chemistry 
comprising  the  compounds  of  metals  with  one  another  has  been  largely  neglected  by 
chemists  in  the  past”  .  For  consistency,  the  element  B  is  chosen  so  that  sBB  is  more 
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negative  than  sAA.  Results  are  shown  in  Figure  5,  and  values  are  tabulated  in 
Supplementary  Infonnation,  Table  S4. 

sAb  values  satisfy  general  expectations.  Since  compounds  are  fonned  in  each  of 
the  binary  systems  studied  here,  a  negative  deviation  from  ideal  solid  solutions  is 
expected.  Consistent  with  this,  sAB  never  exceeds  (sAA  +  sBB)/ 2,  the  criterion  for  ideal 

A  B  A  B 

solid  solutions.  sAB  also  does  not  exceed  the  weighted  average  of  sM  y  and  sBB  v  (Al- 
Nb  is  an  exception)  and  it  intersects  this  weighted  average  at  the  atom  fraction  of 
element  B,  fB,  of  the  most  solute-rich  compound  in  that  binary  system.  sAB  usually  falls 
between  sAABy  and  £BBBy ,  but  is  sometimes  more  negative  than  both  of  these  values. 
Finally,  sAB  values  calculated  from  Equation  9  using  the  ‘elemental’  standard  state 
(commonly  used  in  classical  thermodynamics,  where  su  =  0)  are  essentially  independent 
of  fB.  We  validate  the  unlike  atom  bond  energies  by  predicting  enthalpies  of  formation 

A  B 

for  ternary  compounds,  AjH(AxByCz),  from  eA  y  and  %  and  from  the  numbers  of  i-j 

ABC 

bonds  per  mole  of  AxByCz,  R  '  y  (see  Methods).  Agreement  is  good,  and  is  generally 
within  experimental  error  (Figure  6).  Additional  data  supporting  this  comparison  is 
given  in  Supplementary  Information,  Table  S5. 

To  our  knowledge,  there  are  very  few  published  values  against  which  these 

2 1 

condensed  bond  enthalpies  can  be  compared.  Atomic  potentials  for  the  Ni-Zr  system  “ 
give  sNiNi  and  S/r/r  that  are  between  one  half  and  two  thirds  of  the  values  determined 
here,  and  give  sNizr  values  that  range  from  0.62  to  1.03  of  the  present  values.  Other 
potentials  seem  to  give  values  that  are  much  lower  .  The  crystal  orbital  Hamiltonian 
population  (COHP),  although  not  a  bond  strength  in  the  manner  discussed  here,  is 
nevertheless  an  indication  of  the  covalent  contribution  to  interatomic  bonding  .  COHP 
values  for  Fe-Fe  bonding  range  from  one  half  to  twice  the  value  reported  here  for  SpeFe, 
and  values  for  metal-metalloid  bonding  range  from  about  -2  eV  to  -3.5  eV  ’  .  Bond 
dissociation  enthalpies  for  diatomic  gaseous  elements  generally  range  from  about  3 
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times  smaller  to  4  times  larger  than  the  Sa  values  for  condensed  elements,  consistent 
with  the  expectation  that  these  values  may  not  be  comparable  5. 

Sab  is  relatively  insensitive  to  composition  when  saa~£bb,  but  shows  a  mild 
dependence  on  fB  when  sBb  is  significantly  more  negative  than  sAA.  In  these  systems,  sAB 
is  typically  more  negative  for  A-rich  compounds  and  becomes  slightly  less  negative 
with  increasing  fB.  It  is  possible  that  B-B  bonds  are  increasingly  preferred  as  sBB 
becomes  more  negative  than  sAA,  but  an  analysis  of  the  fractions  of  A-A,  A-B  and  B-B 
bonds  in  the  studied  intennetallic  structures  does  not  conclusively  support  this  idea.  The 
fractions  of  bonds  between  like  and  unlike  atom  pairs  are,  however,  consistently 
influenced  by  atom  size.  Relative  to  the  fractional  number  of  bonds  expected  for  an 
ideal,  random  binary  solution  of  equal-sized  atoms  9,  there  are  fewer  bonds  between  the 
smaller  atoms  across  the  entire  range  of  fB,  there  are  more  bonds  between  unlike  atoms 
for  compounds  rich  in  the  smaller  atom,  and  there  are  more  bonds  between  the  larger 
atoms  for  compounds  rich  in  the  larger  atom  (see  Supplementary  Information, 

Figure  SI).  These  trends  are  consistent  with  simple  topological  arguments  and  with  a 
general  size  dependence  acknowledged  in  conventional  thermodynamics.  The  present 
work  gives  explicit  relations  between  numbers  of  bonds  and  £y,  but,  to  our  knowledge, 
the  relationships  between  atomic  structure  (and  hence  the  frequency  of  A-A,  A-B  and 
B-B  bonds),  topology  (atom  size)  and  the  relative  values  of  sAA,  sAB  and  sBb  are  not 
clearly  understood  and  bear  further  study. 

The  present  approach  offers  the  simplicity  of  pair  interactions,  but  higher-order 
energy  terms  are  embedded  in  the  bond  enthalpies  reported  here,  since  changes  in  total 
enthalpies  are  used  to  calculate  £y.  The  magnitude  of  higher-order  bond  interactions  can 
be  estimated  through  analysis  of  stacking  fault  energies,  which  arise  from  a  change  in 
second-neighbour  bonding.  Assigning  measured  stacking  fault  energies  to  the  bonds 
across  a  stacking  fault  defect  and  normalizing  by  sti  shows  that  second-neighbour 
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interactions  typically  range  from  0.5-5%  of  e-u  for  fee  metallic  elements  and  from  6-16% 
for  hep  metals  (see  Supplementary  Information,  Table  S6).  Measurement  errors  for 
Af/I(A,gas)  and  AjH(AxBy)  give  an  average  uncertainty  in  sy  of  about  ±3%,  and  this  is 
suggested  as  a  basic  limit  in  precision  of  Sy  values  given  here.  For  unusual  cases  where 
the  first  coordination  shell  of  the  relevant  crystal  structure  is  ambiguous  (see 
Methods),  Sy  errors  can  be  as  large  as  ±10%. 

Standard  treatments  of  solution  models  emphasize  the  difference  between  sab  and 
the  average  of  sAa  and  sBB,  but  the  importance  of  relative  differences  between  eaa  and 
£BB  are  rarely  discussed.  The  ‘elemental’  standard  state  does  not  allow  distinction 
between  eaa  and  sBB,  while  this  difference  can  be  quantified  with  the  gas  standard  state. 
The  structure  and  stability  of  condensed  phases  may  depend  on  the  relative  magnitudes 
of  Eaa,  £ab  and  ebb,  so  that  the  gas  standard  state  used  here  for  Sy  is  likely  to  be  useful 
for  understanding  complex,  multi-component  substances  including  ‘superalloys’  that 
can  contain  as  many  as  a  dozen  elements,  high  entropy  alloys  and  metallic 
glasses  27’28. 

Although  the  current  approach  is  applied  here  only  to  metallic  systems,  it  may  be 
more  generally  useful  for  other  condensed,  inorganic  substances  including  compounds 
of  metal/semi-metal  and  metal/non-metal  atoms.  It  may  also  be  useful  in  condensed 
materials  where  ionic  bonding  dominates,  as  ionic  bonding  is  non-directional  and  often 
produces  efficiently  packed  atomic  structures  with  high  coordination  numbers.  It  is  not 
certain  at  present  if  covalent  materials  may  be  treated  with  appropriate  adaptations  of 
the  methods  developed  here.  Adjustment  of  bond  energies  based  on  the  local  atomic 
coordination  in  the  substance  of  interest  gives  good  results  in  the  present  work  for 
extending  bond  enthalpies  from  elements  to  binary  compounds  and  from  binary  to 
ternary  compounds,  and  is  a  suggested  approach. 
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Condensed  bond  enthalpies  are  vital  companions  to  bond  dissociation  enthalpies, 
and  may  provide  a  similar  impact  to  the  understanding  of  liquid  and  solid  substances 
and  their  reactions.  Condensed  bond  enthalpies  can  give  essential  insights  into  the 
structure  and  stability  of  complex  systems,  such  as  ternary  and  higher  order  inorganic 
compounds,  amorphous  metals,  and  high  entropy  alloys,  where  structural  complexity 
challenges  current  modeling  approaches,  and  where  insights  are  likely  to  come  from 
trends  established  by  the  study  of  many  systems  rather  than  detailed  investigations  in  a 
small  number  of  selected  systems.  Condensed  bond  enthalpies  are  also  expected  to  give 
new  insights  into  the  kinetics  of  chemical  reactions  such  as  catalysis  and  the  fragility  of 
liquids  ~  ,  where  the  relative  magnitudes  of  bond  energies  at  a  local  atomic  scale  are 
important.  Condensed  bond  enthalpies  provide  new  data  to  support  the  established  field 
of  solution  thennodynamics  and  may  give  new  insights  by  quantifying  differences 
between  saa  and  Sbb ■  Finally,  observations  in  the  present  work  link  chemistry  and 
structural  topology,  suggesting  an  intriguing  approach  to  combine  topology  and  energy 
of  atomic  structures  in  liquids  and  solids. 

METHODS 

Enthalpy  of  broken  bonds.  Cohesive  energy  is  intimately  connected  to  the  spatial 
distribution  of  electrons  .  The  redistribution  of  electrons  at  defects  is  thus  expected 
to  alter  the  enthalpy  of  the  broken  bonds  associated  with  those  defects.  Consider  the 
redistribution  of  electrons  associated  with  bond  breaking:  electrons  in  the  initial  bond 
become  partitioned  between  the  broken,  or  ‘dangling’,  bond  and  the  bonds  that  remain 
intact  around  the  defect  site.  We  propose  a  simple  estimate  of  this  partitioning,  where 
the  fraction  of  enthalpy  remaining  in  the  broken  bond  is  equal  to  the  fractional 
coordination  number  of  atoms  in  the  first  coordination  shell  of  the  defect.  This 
correction  is  simplest  for  free  surfaces.  Surface  atoms  on  (1 1  l)fCC  planes  have  a 
coordination  of  9  and  an  initial  coordination  number  of  12,  so  that  the  fractional 
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enthalpy  remaining  in  the  broken  bonds  is  3/4.  The  fraction  is  2/3  for  (100)fCC  planes  and 
3/4  for  (110)  fCC  planes,  giving  an  average  fractional  enthalpy  of  13/18  (about  0.722)  for 
broken  bonds  associated  with  low-index  fee  planes.  Fractional  coordinations  for  surface 
atoms  are  5/7  (about  0.714)  for  each  of  the  (100)bcc,  (1 10)bcc  and  (1 12)bcc  planes. 

The  correction  for  the  enthalpy  of  broken  bonds  at  a  vacancy  requires  terms  for 
atoms  surrounding  both  the  site  where  the  atom  is  removed  and  where  it  is  placed  on  the 
surface.  The  twelve  atoms  surrounding  a  vacant  site  in  an  fee  structure  have  a  final 
coordination  of  1 1  and  an  initial  coordination  of  12.  Placing  the  displaced  atom  in  a 
surface  vacancy  on  a  (1 1  l)fCC  plane  changes  its  coordination  number  from  12  to  9, 
changes  the  coordination  number  for  the  three  atoms  at  the  bottom  of  the  surface 
vacancy  from  11  to  12,  and  changes  the  coordination  number  from  8  to  9  for  the  six 
remaining  atoms  that  form  the  surface  vacancy.  The  average  change  in  coordination 
number  for  the  atoms  associated  with  this  defect  is  thus 

12 (11 /12)+1  (9/12 )+3 (12/11) +6(9/8)  _  Q  ^ 

Similar  analysis  gives  fractions  of  0.981  for  (100)fcc  planes  and  0.970  for  (1 10)fcc  planes, 
for  an  average  of  0.980  for  low-index  surfaces  in  the  fee  structure.  Average  corrections 
are  0.981  for  ( 1 00)bcc,  (11 0)bcc  and  (112)bcc  planes  and  0.981  for  (0001)ftcp,  (1010)hcp, 
(1011) hcp  and  (1121) hcp.  planes. 

A  similar  correction  is  made  to  account  for  the  bonds  broken  upon  melting.  The 
change  in  coordination  occurs  equally,  in  a  stochastic  sense,  for  all  atoms.  Only  fee  and 
hcp  structures  are  considered  here,  so  that  the  initial  coordination  number  is  12,  the  final 
coordination  number  is  Z-AZ  =11.3  and  the  fractional  enthalpy  retained  in  the  broken 
bonds  is  0.942. 
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Temperature  adjustment  of  Saa-  The  enthalpy  content  of  a  solid  at  298  K  in  the 
present  work,  -A ,fH(A,gas),  becomes  less  negative  as  the  solid  is  heated  above  this 
temperature.  This  enthalpy  reduction  is  (. Hjm-H298 )  at  the  melting  temperature,  so  that 
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£aa  =  £AA 


298) 


&fH(A,gas  )  J 


(M2) 


1 8 

When  available,  tabulated  values  of  (HT  -H29s)  are  used  in  the  present  analyses.  When 


they  are  not  available,  sAA  is  estimated  as 


£aa  =  £aa 


1  - 


Cp  (Tm  —298) 


AfH(A,gas)j 


(M3) 


where  Cp  is  the  elemental  heat  capacity. 

Surface  energies.  The  surface  energy,  y,  is  given  as 


y  -  Fs  -  Es  -  TSs 


(M4) 


where  F,  E  and  S  are  the  Helmholtz  free  energy,  internal  energy  and  entropy  and  the 
subscript,  s,  indicates  properties  at  a  free  surface.  The  internal  energy  is  estimated  as 


Es  =  -E(e^T  )  =  — E (£aa) 


1  - 


( HTn 


_#298) 


AfH(A,gas ) 


(M5) 


T 

where  E  is  the  number  of  bonds  per  unit  area.  £AAeas  is  the  bond  enthalpy  at  the 
temperature  which  the  surface  energy  is  measured,  Tmeas,  and  is  detennined  as  described 
above  for  eAA  .  The  surface  entropy  is  evaluated  by  taking  the  temperature  derivative  of 
Fs  and  rearranging  terms 


(E)(£2L4  ) 

d(#r~#298) 

_dt± 

AfH{A,gas ) 

dT 

dT 

(M6) 


The  temperature  dependence  of  ( HT  —  H2 9s)  is  detennined  using  tabulated  data  18  from 
temperatures  just  above  (T2)  and  just  below  (Tl)  the  measurement  temperature,  so  that 
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(E)Oa4  ) 

(tf  7-2  — 298  )— (W  7-1  298  )] 

1 

a. 

A  fH(A,gas) 

T2—T1  \ 

dT 

(M7) 


Measured  values  of  d  Fs/dT  are  given  in  15.  The  surface  energy  is  calculated  by 
inserting  Equations  M5  and  M7  into  Equation  M4.  The  value  of  saa  used  here  includes 
correction  for  the  change  in  coordination  number  associated  with  bonds  broken  upon 
melting,  as  described  earlier  in  Methods.  Tabulated  values  of  E  (Table  S2)  and  of  data 
used  in  the  calculation  of  y  (Table  S3)  are  given  in  Supplementary  Information. 


Calculation  of  Sab-  Bond  counting  and  adjustment  of  ea  to  account  for  the  different 
structures  in  the  element  and  the  compound  are  used  to  transform  Equation  8  for  sAB 
into  Equation  9.  Crystallographic  data  enable  bond  counting.  The  Pearson  symbol  for 
the  equilibrium  structure  ’  points  to  detailed  information  of  the  unit  cell  .  The 
numbers  of  A-A,  A-B  and  B-B  bonds  between  the  central  atom  and  each  of  the  atoms 
in  the  first  coordination  shell  are  counted  for  each  of  the  constituent  polyhedra  in  the 
unit  cell.  These  values  are  multiplied  by  appropriate  site  occupancies  and  site 
multiplicities,  and  the  sums  are  halved  to  avoid  double-counting  of  bonds.  The  resulting 

A  B 

numbers  of  bonds  per  unit  cell  of  AXBV,  pA  y  are  converted  to  the  number  bonds  per 
mole  of  AxBy,  P^xBy  via 

(M8) 

Here  x  and  v  are  the  numbers  of  A  and  B  atoms,  respectively,  in  the  AxBy  fonnula  unit 
and  U  is  the  number  of  atoms  per  unit  cell.  While  polyhedral  coordination  numbers  in  34 
are  usually  reliable,  in  rare  cases  apparent  discrepancies  exist.  We  perform  a  critical 
assessment  of  the  coordination  number  for  each  polyhedron  to  include  atoms  up  to  a 
distance  of  1 .25  times  the  minimum  atomic  separation  for  each  of  the  A-A,  A-B  and 
B-B  atom  pairs.  Where  a  discrepancy  between  reported  and  assessed  coordination 
numbers  is  found,  we  use  the  average  of  the  two  values  with  an  associated  error  in 

A  B 

pA  y .  This  results  in  a  larger  error  for  sab- 


pAxBy 

ij 


_  pH 
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We  propose  that  bond  enthalpy  depends  on  structure  via  the  number  of  bonds 
formed  per  atom,  as  suggested  by  Figure  lb.  Thus,  the  more  bonds  fonned  per  atom,  the 
lower  is  the  enthalpy  per  bond.  We  apply  a  simple  estimate  of  this  effect 


,AxBy 


(M9) 


—elem 


where  pl“  —  %aa  /2  is  the  average  number  of  bonds  per  atom  in  the  pure  element,  and 
—AxBy  _  is  the  average  number  of  bonds  per  atom  in  the 


— elem 


—ArBv 


AXBV  structure,  p  is  usually  either  6  (for  fee  and  hep)  or  7  (for  bee),  while  p  *  y  is 


typically  between  6  and  7.  A  similar  correction  is  applied  for  ternary  compounds, 


AXBy  Cz 
u 


AXBy  Cz 
F 

CAB 


AXBy  Cz 
F 

cAC 


(MlOa) 


—AXCZ 


—AXByCZ 


) 


£ac> 


AXBy  Cz 
F 

cBC 


(=W7 )^c  (Ml Ob) 


_ A  ft 

where  p  *  y  is  determined  for  the  binary  compound  at  the  same  ratio  of  x:y  as  in  the 
ternary,  or  as  the  weighted  average  of  the  binary  compounds  that  bracket  this 
composition,  sab  is  determined  by  inserting  fB  =  x/(x+y)  into  a  linear  regression  of  sab 
vs.  Jb  from  the  appropriate  binary  systems. 


Enthalpies  of  ternary  compounds.  The  formation  enthalpy  of  a  ternary  compound, 

AxByCz,  is 


A fH(AxByCz')  —  H  {AxByCz)  —  xH(A,  cond)  —  yH{B,cond)  —  zH(C,cond)  (Mil) 
Expanding  H{AxByCz)  and  substituting  —Af  H(i,  gas)  for  H (i,  cond )  gives 


AfH{AxByCz)  — 


AXByCZ  AXByCZ 

"aa  £aa 

AXByCZ  AXByCZ 
, BC  eBC 


Ax  By  Cz  Ax  By  Cz 

+  "ab  £ab 

Ax  By  Cz  Ax  By  Cz 

+  "cc  £cc 


AXByCZ  AXByCZ  AXByCZ  AXByCZ 
+  *AC  £AC  +  ^BB  £BB  + 

+  xA fH{A,  gas )  +  yAfH(B,  gas )  +  zAfH{C,  gas ) 


(Ml  2) 
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Table  1  Elemental  condensed  bond  enthalpies,  sAA 


Elemen 

Atomic 

Numbei 

<  > 

$  S, 

Elemen 

Atomic 

Numbei 

<  > 

$  S. 

Elemen 

Atomic 

Numbei 

<  > 

$  S. 

H 

1 

-2.26±0.00001 

As 

33 

-1.05±0.045 

Tb 

65 

-0.671±0.004 

Li 

3 

-0.236±0.002 

Se 

34 

-1.57±0.028 

Dy 

66 

-0.502±0.004 

Be 

4 

-0.560±0.009 

Br 

35 

-1 . 1 6±0.00 1 

Ho 

67 

-0.51 9±0. 004 

B 

5 

-1.80±0.016 

Rb 

37 

-0.120±0.001 

Er 

68 

-0.547±0.004 

C 

6 

-4.95±0.003 

Sr 

38 

-0.283±0.003 

Tm 

69 

-0.401±0.004 

N 

7 

-4.90±0.004 

Y 

39 

-0.734±0.004 

Yb 

70 

-0.269±0.004 

0 

8 

-2.58±0.001 

Zr 

40 

-1.05±0.015 

Lu 

71 

-0.739±0.004 

F 

9 

-0.823±0.003 

Nb 

41 

-1.09±0.012 

Hf 

72 

-1.07±0.011 

Na 

11 

-0.159±0.001 

Mo 

42 

-0.976±0.006 

Ta 

73 

-1.16±0.004 

Mg 

12 

-0.254±0.001 

To 

43 

-1.17±0.023 

W 

74 

-1.26±0.009 

Al 

13 

-0.572±0.007 

Ru 

44 

-1.12±0.011 

Re 

75 

-1.34±0.011 

Si 

14 

-2.33±0.042 

Rh 

45 

-0.960±0.007 

Os 

76 

-1.36±0.011 

P 

15 

-2.19±0.007 

Pd 

46 

-0.651±0.004 

Ir 

77 

-1.16±0.007 

S 

16 

-2.87±0.002 

Ag 

47 

-0.492±0.001 

Pt 

78 

-0.977±0.002 

Cl 

17 

-1.26±0.0001 

Cd 

48 

-0.193±0.003 

Au 

79 

-0.636±0.004 

K 

19 

-0.132±0.001 

In 

49 

-0.420±0.007 

Hg 

80 

-0.106±0.0001 

Ca 

20 

-0.307±0.001 

Sn 

50 

-0.624±0.003 

TI 

81 

-0.315±0.001 

Sc 

21 

-0.653±0.013 

Sb 

51 

-0.914±0.009 

Pb 

82 

-0.337±0.001 

Ti 

22 

-0.817±0.005 

Te 

52 

-1.36±0.015 

Bi 

83 

-0.483±0.005 

V 

23 

-0.763±0.012 

1 

53 

-1.11  ±0.004 

Po 

84 

-0.491±0.010 

Cr 

24 

-0.589±0.006 

Cs 

55 

-0.1 13±0.002 

Ra 

88 

-0.235±0.005 

Mn 

25 

-0.448±0.007 

Ba 

56 

-0.265±0.007 

Ac 

89 

-0.701±0.014 

Fe 

26 

-0.61 5±0. 002 

La 

57 

-0.745±0.004 

Th 

90 

-1.04±0.010 

Co 

27 

-0.737±0.015 

Ce 

58 

-0.726±0.004 

Pa 

91 

-0.834±0.017 

Ni 

28 

-0.743±0.015 

Pr 

59 

-0.617±0.004 

U 

92 

-0.921±0.014 

Cu 

29 

-0.583±0.002 

Nd 

60 

-0.565±0.004 

Np 

93 

-0.642±0.013 

Zn 

30 

-0.225±0.001 

Sm 

62 

-0.357±0.004 

Pu 

94 

-0.51 1  ±0.01 0 

Ga 

31 

-0.805±0.006 

Eu 

63 

-0.263±0.003 

Ge 

32 

-1.93±0.016 

Gd 

64 

-0.687±0.004 
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Figure  1  |  Elemental  condensed  bond  enthalpies  vs  melting  temperature,  Tm.  a 
Condensed  bond  enthalpies  calculated  using  coordination  numbers  of  the 
equilibrium  structures,  b  Condensed  bond  enthalpies  calculated  using  a  space¬ 
filling  coordination  number  of  12  for  all  elements. 
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Figure  2  |  Measured  and  predicted  monovacancy  formation  energies  in  metallic 
elements.  The  solid  line  is  a  least-squares  fit  and  the  dotted  line  represents 
perfect  agreement.  The  error  bars  give  the  reported  range  in  experimental 
values  for  the  elements  indicated.  Bond  enthalpy  estimates  exceed  measured 
values  by  about  0.1  eV,  within  typical  experimental  error.  Experimental  data  are 
from  35'36’37'38. 
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Figure  3  |  Measured  and  predicted  enthalpies  of  fusion  for  metallic  elements. 
The  dotted  line  represents  perfect  agreement.  Error  bars  show  ranges  in 
experimental  values,  taken  from  4'5. 
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Figure  4  |  Comparison  of  measured  and  predicted  solid  surface  energies,  y. 
Error  bars  for  predicted  data  give  the  range  in  values  predicted  for  the  different 
crystallographic  planes.  A  modest  error  of  ±5%  is  assumed  for  measured 
values.  The  solid  line  is  a  linear  regression  for  elements  where  measurements 
are  made  at  <  1723  K,  and  the  dotted  line  represents  perfect  agreement. 
Measurements  for  Nb,  Mo  and  W  are  made  at  much  higher  temperatures, 
where  impurity  adsorption  may  alter  surface  energies.  Experimental  data  are 
taken  from  15. 
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Figure  5  |  For  each  binary  system  are  plotted:  saa  for  element  A  (at  atom  fraction  B  =  0)  and  in  the 
binary  compounds;  sBb  for  element  B  (at  atom  fraction  B  =  1 )  and  in  the  binary  compounds;  and 
sab  relative  to  both  the  gas  and  ‘elemental’  standard  states.  The  compositionally-weighted 

A  B  A  B  _ 

averages  of  ej  y  and  eBB  y  are  shown  by  open  crosses.  The  more  negative  sn  of  the  two 
elements  is  B  in  all  plots. 
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Figure  6  |  Comparison  of  measured  and  predicted  enthalpies  of  formation  for 

A  B 

ternary  compounds,  AfH(AxByCz).  Predictions  are  made  using  eax  y  and  with  £,j 
values  relative  to  both  the  gas  and  metal  standard  states.  The  dotted  line 
represents  perfect  agreement.  Experimental  data  are  taken  from  39. 
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SUPPLEMENTARY  INFORMATION 

Table  SI  Elemental  bond  enthalpies,  sAa 


Element 

* 

't/T  L_ 
re  o 

3  E 

X  “s 
<  3 

Q_ 

<  > 

5  3 

Element 

* 

'c/T  Ti- 
ra  o 

3  E 

X  “5 

<  3 

2 

Q_ 

5  > 

<  <D 

05 

H 

217.998±0.0006 

i 

-2.26±0. 00001 

Cd 

1 1 1 ,8±0.2 

6 

-0.193±0.003 

Li 

159. 3±1 .0 

7 

-0.236±0.002 

In 

243±4 

6 

-0.420±0.007 

Be 

324±5 

6 

-0.560±0.009 

Sn 

301 ,2±1 .5 

5 

-0.624±0.003 

B 

565±5 

3.25 

-1.80±0.016 

Sb 

264.4±2.5 

3 

-0.914±0.009 

C 

716.68±0.45 

1.5 

-4.95±0.003 

Te 

196.6±2.1 

1.5 

-1.36±0.015 

N 

472.68±0.4 

1 

-4.90±0.004 

1 

106.76±0.04 

1 

-1.11  ±0.004 

O 

249. 229±0. 0002 

1 

-2.58±0.001 

Cs 

76.5±1 .0 

7 

-0.1 13±0.002 

F 

79.38±0.3 

1 

-0.823±0.003 

Ba 

179.1±5.0 

7 

-0.265±0.007 

Na 

107.5±0.7 

7 

-0.159±0.001 

La 

431  ±2.1 

6 

-0.745±0.004 

Mg 

147.1±0.8 

6 

-0.254±0.001 

Ce 

420.1±2.1 

6 

-0.726±0.004 

Al 

330. 9±4 

6 

-0.572±0.007 

Pr 

356.9±2.1 

6 

-0.61 7±0. 004 

Si 

450±8 

2 

-2.33±0.042 

Nd 

326.9±2.1 

6 

-0.565±0.004 

P 

316.5±1.0 

1.50 

-2.19±0.007 

Sm 

206.7±2.1 

6 

-0.357±0.004 

S 

277.17±0.15 

1 

-2.87±0.002 

Eu 

177.4±2.1 

7 

-0.263±0.003 

Cl 

121 ,301±0.008 

1 

-1.26±0.0001 

Gd 

397.5±2.1 

6 

-0.687±0.004 

K 

89±0.8 

7 

-0.132±0.001 

Tb 

388.7±2.1 

6 

-0.671±0.004 

Ca 

177.8±0.8 

6 

-0.307±0.001 

Dy 

290.4±2.1 

6 

-0.502±0.004 

Sc 

377. 8±4 

6 

-0.653±0.013 

Ho 

300.6±2.1 

6 

-0.51 9±0. 004 

Ti 

473±3 

6 

-0.81 7±0. 005 

Er 

316.4±2.1 

6 

-0.547±0.004 

V 

515. 5±8 

7 

-0.763±0.012 

Tm 

232.2±2.1 

6 

-0.401±0.004 

Cr 

397.48±4.2 

7 

-0.589±0.006 

Yb 

155.6±2.1 

6 

-0.269±0.004 

Mn 

283.3±4.2 

6.55 

-0.448±0.007 

Lu 

427.6±2.1 

6 

-0.739±0.004 

Fe 

415.5±1.3 

7 

-0.61 5±0. 002 

Hf 

618.4±6.3 

6 

-1.07±0.01 1 

Co 

426.7±8.5 

6 

-0.737±0.015 

Ta 

782±2.5 

7 

-1.16±0.004 

Ni 

430.1±8.4 

6 

-0.743±0.015 

W 

851  ±6.3 

7 

-1.26±0.009 

Cu 

337. 4±1 .2 

6 

-0.583±0.002 

Re 

774±6.3 

6 

-1.34±0.01 1 

Zn 

130.4±0.4 

6 

-0.225±0.001 

Os 

787±6.3 

6 

-1.36±0.01 1 

Ga 

271 ,96±2.1 

3.5 

-0.805±0.006 

Ir 

669±4 

6 

-1.16±0.007 

Ge 

372±3 

2 

-1.93±0.016 

Pt 

565. 7±1 .3 

6 

-0.977±0.002 

As 

302.5±13 

3 

-1.05±0.045 

Au 

368.2±2.1 

6 

-0.636±0.004 

Se 

227. 2±4 

1.5 

-1.57±0.028 

Hg 

61 ,38±0.04 

6 

-0.106±0.0001 

Br 

1 1 1 ,87±0.1 

1 

-1.16±0.001 

TI 

182.2±0.4 

6 

-0.31 5±0. 001 

Rb 

80.9±0.8 

7 

-0.120±0.001 

Pb 

195.2±0.8 

6 

-0.337±0.001 

Sr 

164±1.7 

6 

-0.283±0.003 

Bi 

209.6±2.1 

4.5 

-0.483±0.005 

Y 

427.4±2.1 

6 

-0.734±0.004 

Po 

142±2.8 

3 

-0.491±0.010 

Zr 

610±8.4 

6 

-1.05±0.015 

Ra 

159±3.2 

7 

-0.235±0.005 

Nb 

733±8 

7 

-1.09±0.012 

Ac 

406±8.1 

6 

-0.701±0.014 

Mo 

658.98±3.8 

7 

-0.976±0.006 

Th 

602±6 

6 

-1.04±0.010 

Tc 

678±13.6 

6 

-1.17±0.023 

Pa 

563±1 1 .3 

7 

-0.834±0.017 

Ru 

650.6±6.3 

6 

-1.1 2±0.0 1 1 

U 

533±8 

6 

-0.921±0.014 

Rh 

556±4 

6 

-0.960±0.007 

Np 

464.8±9.3 

7.5 

-0.642±0.013 

Pd  376.6±2.1  6  -0.651±0.004 

Ag  284.9±0.8  6  -0.492±0.001 

'  Assessed  from  values  reported  in  4'5'11'12. 

Pu 

345±6.9 

7 

-0.51 1  ±0.01 0 
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Table  S2  Bonds  broken  per  unit  area,  E 


Plane  and 

structure 

E 

(100)fcc 

1 

2r2 

(110)fcc 

3 

4V2r2 

(lll)fcc 

3 

4V3 r2 

(1  00)bcc 

9 

16r2 

(110)bCc 

3 

4V2 r2 

(112)bCc 

V3 

2V2r2 
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Table  S3  Measured  and  calculated  surface  energies 


Element 

a 

Ymeas 

(mJ  nr2) 

^  £ 

(/>  “> 

%  £ 

* 

if) 

ro 

0) 

E 

1- 

OO  1 
ci  — 

CM  O 

^  E 
£3 

H 

<1 

OO* 

O  V 

£  * 

X  1 

<  3 

Yioo 

(mJ  m-2) 

Yiio 

(mJm5) 

Ym 

(mJm5) 

Y112 

(mJm5) 

Ag 

1100 

-0.47 

1223 

26077 

31.38 

833 

1071 

769 

Au 

1400 

-0.43 

1273 

27041 

30.54 

1339 

1645 

1243 

Cu 

1780 

-0.50 

1198 

24606 

29.50 

1759 

2117 

1618 

Fe 

1950 

-0.90 

1723 

55094 

41.42 

1641 

1458 

1924 

Co 

1970 

-0.52  b 

1627 

49064 

39.75 

2329 

2764 

2101 

Nb 

2100 

-0.52 b 

2523 

65380 

33.89 

3150 

2894 

3546 

Pt 

2200 

-0.60 

1573 

37319 

32.64 

2219 

2711 

2040 

Ni 

2280 

-0.55 

1333 

33195 

35.35 

2286 

2739 

2102 

Cr 

2300 

-0.52 b 

1673 

44590 

42.26 

2128 

1956 

2394 

Mo 

2400 

-0.20 

2623 

73976 

40.17 

3759 

3514 

4139 

W 

2800 

-0.52 b 

2273 

55679 

31.38 

3888 

3598 

4337 

a  From  15). 

b  Estimated  as  the  average  of  the  measured  values. 
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Table  S4  Calculated  values  of  8Ab 


Compound 

AXBy 

Pearson 

Symbol 

(Prototype) 

A-A  bonds 
per  unit  cel 

A-B  bonds 
per  unit  cel 

B-B  bonds 
per  unit  cel 

CQ  ® 

Jf  0 

i  E 

X  -3 

<  i 

cT  .0 

^  3  > 

a 

T3 

oT  .0 

*  “a  > 

^  a 

■0 

c 

0 

±2 

m 

£  a 

AI3Hf 

til 6  (AI3Zr) 

48 

48 

4 

-168±14 

-0.549±0.021 

-1.03±0.039 

-0.902±0.090 

AI2Hf 

hP12  (MgZn2) 

24 

48 

8 

-1 44±1 2 

-0.514±0.007 

-0.961±0.01 1 

-0.813±0.018 

AI3Hf2 

oF40  (AI3Zr2) 

48 

148 

48 

-240±20 

-0.562±0.026 

-1.05±0.049 

-0.860±0.081 

AI3Hf4 

hP7  (AI3Zr4) 

6 

24 

17 

-322±27 

-0.511  ±0.007 

-0.955±0.01 1 

-0.832±0.019 

AI2Hf3 

tP20  (AI2Zr3) 

8 

64 

59 

-220±18 

-0.524±0.019 

-0.979±0.036 

-0.805±0.070 

AIHf2 

tl  1 2  (AI2Cu) 

4 

32 

44 

-1 23±1 0 

-0.514±0.007 

-0.961±0.01 1 

-0.804±0.019 

AIi2Mo 

cl26  (Ah2W) 

12 

1 

13 

-195±87 

-0.619±0.007 

-1.23±0.006 

-1.07±0.085 

AI5Mo 

hP12  (AI5W) 

5 

1 

6 

-192±35 

-0.572±0.007 

-1.14±0.006 

-1.02±0.037 

AI4M0 

mC30  (AI4W) 

4 

1 

5 

-185±29 

-0.544±0.025 

-1.08±0.049 

-1.03±0.092 

AIsMo3 

mC22  (AI8Mo3) 

8 

3 

11 

-484±41 

-0.576±0.007 

-1.15±0.009 

-1.02±0.036 

AIMo3 

cP8  (Cr3Si) 

1 

3 

4 

-142±15 

-0.508±0.007 

-1.01±0.006 

-0.851±0.019 

AI3Nb 

tl8  (AI3Ti) 

24 

24 

0 

-132±1 1 

-0.572±0.007 

-1.27±0.012 

-1.03±0.020 

AINb2 

tP30  (CrFe) 

12 

94 

88 

-75±6 

-0.530±0.021 

-1.18±0.047 

-0.896±0.077 

AINb3 

cP8  (Cr3Si) 

0 

24 

30 

-76±6 

-0.508±0.007 

-1.13±0.012 

-0.844±0.015 

AI3Ni 

0PI6  (CFe3) 

68 

34 

2 

-192±16 

-0.573±0.020 

-0.686±0.025 

-0.874±0.096 

AI3Ni2 

hP5  (AI3Ni2) 

9 

16 

3 

-306±26 

-0.612±0.007 

-0.796±0.015 

-0.905±0.029 

AINi 

cP2  (CsCI) 

3 

8 

3 

-132±1 1 

-0.490±0.007 

-0.637±0.015 

-0.734±0.022 

AI3Ni5 

oCI  6  (Ga3Pt5) 

8 

56 

32 

-435±37 

-0.572±0.007 

-0.743±0.015 

-0.81 8±0. 024 

AINi3 

cP4  (AuCu3) 

0 

12 

12 

-1 64±1 4 

-0.572±0.007 

-0.743±0.015 

-0.799±0.023 

AI3Pt2 

hP5  (AI3Ni2) 

9 

16 

3 

-435±37 

-0.612±0.007 

-1.05±0.022 

-1.12±0.029 

AlPt 

cP8  (FeSi) 

12 

28 

12 

-204±17 

-0.528±0.007 

-0.902±0.022 

-1.02±0.029 

AI3Pt5 

0PI6  (Ge3Rh5) 

5 

55 

38 

-704±59 

-0.560±0.042 

-0.957±0.072 

-0.993±0.081 

AIPt3 

tP16  (GaPt3) 

0 

48 

50 

-280±24 

-0.560±0.007 

-0.957±0.022 

-0.996±0.025 

AI3Ti 

tl8  (AI3Ti) 

24 

24 

0 

-147±3 

-0.572±0.007 

-0.817±0.005 

-0.821±0.008 

AITi 

tP4  (AuCu) 

4 

16 

4 

-74±2 

-0.572±0.007 

-0.817±0.005 

-0.790±0.009 

AITi3 

hP8  (Ni3Sn) 

0 

24 

24 

-100±8 

-0.572±0.007 

-0.817±0.005 

-0.781±0.013 

AI3Y 

hP8  (Ni3Sn) 

24 

24 

0 

-190±16 

-0.572±0.007 

-0.734±0.004 

-0.817±0.019 

ai2y 

cF24 (Cu2Mg) 

48 

96 

16 

-1 51  ±1 3 

-0.514±0.007 

-0.660±0.004 

-0.702±0.016 

AIY 

0C8  (BCr) 

4 

28 

20 

-90±8 

-0.528±0.007 

-0.677±0.004 

-0.693±0.017 

ai2y3 

tP20  (AI2Zr3) 

8 

64 

59 

-200±17 

-0.524±0.019 

-0.672±0.025 

-0.699±0.048 

aiy2 

oP12  (Co2Si) 

0 

40 

32 

-105±9 

-0.572±0.007 

-0.734±0.004 

-0.745±0.015 

AI3Zr 

til 6  (AI3Zr) 

48 

48 

4 

-1 63±1 4 

-0.549±0.021 

-1.01±0.039 

-0.892±0.089 

AI2Zr 

hP12  (MgZn2) 

24 

48 

8 

-1 37±1 2 

-0.514±0.007 

-0.948±0.015 

-0.802±0.018 

AI3Zr2 

oF40  (AI3Zr2) 

48 

148 

48 

-235±20 

-0.562±0.026 

-1.04±0.049 

-0.853±0.082 

AlZr 

0C8  (BCr) 

4 

28 

20 

-89±7 

-0.528±0.007 

-0.973±0.015 

-0.755±0.018 

AI2Zr3 

tP20  (AI2Zr3) 

8 

64 

59 

-192±16 

-0.524±0.019 

-0.965±0.036 

-0.783±0.072 

AIZr2 

hP6  (lnNi2) 

0 

22 

14 

-100±8 

-0.572±0.007 

-1.05±0.015 

-0.885±0.018 
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AIZr3 

cP4  (AuCu3) 

0 

12 

12 

-108±9 

-0.572±0.007 

-1.05±0.015 

-0.906±0.019 

CueLa 

oP28  (CeCue) 

106 

76 

0 

-79±1 1 

-0.538±0.017 

-0.687±0.022 

-0.632±0.030 

CusLa 

hP6  (CaCu5) 

21 

18 

0 

-75±3 

-0.538±0.021 

-0.687±0.036 

-0.635±0.004 

Cu2La 

hP3  (AIB2) 

3 

12 

4 

-51  ±3 

-0.552±0.021 

-0.705±0.036 

-0.626±0.005 

CuLa 

oP8  (BFe) 

4 

28 

20 

-32±3 

-0.538±0.021 

-0.687±0.036 

-0.617±0.006 

HfNis 

cF24  (AuBe5) 

96 

64 

0 

-252±21 

-0.669±0.015 

-0.961±0.01 1 

-0.954±0.023 

Hf2Ni7 

mC36  (Ni7Zr2) 

112 

112 

14 

-468±39 

-0.674±0.017 

-0.970±0.024 

-0.950±0.073 

a-HfNia 

hR12  (BaPb3) 

36 

36 

4.5 

-228±19 

-0.699±0.039 

-1.01±0.056 

-1.02±0.124 

Hf7Niio 

oC68  (NiioZr7) 

80 

258 

88 

-1071  ±90 

-0.712±0.035 

-1.02±0.050 

-0.989±0.088 

HfNi 

oC8  (BCr) 

4 

28 

20 

-130±1 1 

-0.686±0.015 

-0.986±0.01 1 

-0.943±0.028 

Hf2Ni 

tl  1 2  (AI2Cu) 

4 

32 

44 

-1 41  ±1 2 

-0.669±0.015 

-0.961±0.01 1 

-0.937±0.026 

LaNis 

hP6  (CaCu5) 

21 

18 

0 

-148±21 

-0.686±0.015 

-0.687±0.004 

-0.772±0.021 

La2Ni7 

hP36  (Ce2Ni7) 

108 

120 

12 

-234±20 

-0.669±0.015 

-0.670±0.004 

-0.750±0.015 

LaNi3 

hR12  (Be3Nb) 

33 

42 

5 

-98±14 

-0.669±0.015 

-0.670±0.04 

-0.742±0.019 

LaNi2 

cF24 (Cu2Mg) 

48 

96 

16 

-75±15 

-0.669±0.015 

-0.670±0.004 

-0.734±0.021 

LaNi 

oC8  (BCr) 

4 

28 

20 

-40±3 

-0.686±0.015 

-0.687±0.004 

-0.745±0.016 

LaaNi 

oP16  (CFe3) 

2 

34 

68 

-52±4 

-0.686±0.025 

-0.687±0.026 

-0.750±0.056 

NbNi3 

oP8  (Cu3Ti) 

24 

24 

0 

-133±1 1 

-0.743±0.015 

-1.27±0.012 

-1.12±0.024 

Nb7Ni6 

hR13  (Fe7W6) 

18 

48 

21 

-281  ±24 

-0.666±0.015 

-1 . 1 4±0.01 2 

-0.979±0.017 

Ni3Ti 

hP16  (Ni3Ti) 

48 

48 

0 

-1 40±1 2 

-0.743±0.015 

-0.817±0.005 

-0.901±0.020 

NiTi2 

cF96  (NiTi2) 

48 

288 

288 

-83±2 

-0.686±0.015 

-0.754±0.005 

-0.81 1±0.01 1 

NiY3 

oP16  (CFe3) 

68 

34 

2 

-76±6 

-0.677±0.025 

-0.686±0.025 

-0.776±0.060 

Ni2Y3 

tP80  (Ni2Y3) 

256 

228 

24 

-1 50±1 3 

-0.693±0.005 

-0.702±0.015 

-0.809±0.033 

NiY 

oP8  (BFe) 

20 

28 

4 

-74±6 

-0.677±0.004 

-0.686±0.015 

-0.794±0.020 

Ni2Y 

cF24 (Cu2Mg) 

16 

96 

48 

-1 1 7±1 0 

-0.660±0.004 

-0.669±0.015 

-0.767±0.017 

Ni3Y 

hR12  (Be3Nb) 

5 

42 

33 

-1 48±1 2 

-0.660±0.004 

-0.669±0.015 

-0.775±0.018 

Ni7Y2 

hR18  (Co7Er2) 

6 

60 

54 

-315±26 

-0.660±0.004 

-0.669±0.015 

-0.775±0.017 

Ni5Y 

hP6  (CaCu5) 

0 

18 

21 

-195±9 

-0.677±0.004 

-0.686±0.015 

-0.795±0.014 

N  i  1 7  Y  2 

hP38  (Nii7Th2) 

0 

76 

164 

-361±30 

-0.697±0.033 

-0.706±0.034 

-0.801±0.087 

Ni5Zr 

cF24  (AuBe5) 

96 

64 

0 

-21 0±1 8 

-0.669±0.015 

-0.948±0.015 

-0.921±0.022 

Ni7Zr2 

mC36  (Ni7Zr2) 

112 

112 

14 

-414±35 

-0.674±0.017 

-0.956±0.024 

-0.926±0.073 

NiioZr7 

oC68  (NiioZr7) 

80 

258 

88 

-884±74 

-0.712±0.035 

-1.01±0.050 

-0.955±0.085 

NiZr 

oC8  (BCr) 

4 

28 

20 

-98±8 

-0.686±0.015 

-0.973±0.015 

-0.892±0.025 

NiZr2 

tl  1 2  (AI2Cu) 

4 

32 

44 

-111  ±9 

-0.669±0.015 

-0.948±0.015 

-0.894±0.023 

b  Assessed  from  values  reported  in  16,17,18,19. 
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Table  S5  Measured  and  calculated  ternary  compound  heats  of  formation 


0) 

0 

0 

0 

o 

O 

O 

o 

N 

O 

>4-J 

•4-* 

'c 

'c 

'c 

'c 

— 

3 

3 

3 

3 

DO 

O 

4_ 

i_ 

4- 

i_ 

X 

-Q 

d) 

0 

0 

0 

< 

E 

Q. 

a 

a 

a 

T3 

CO 

CO 

CO 

CO 

C 

~o 

~o 

~o 

~o 

3 

c  >, 

c 

c 

c 

c 

O 

o 

o 

o 

o 

o 

Q. 

CO  o 

-Q 

-Q 

-Q 

_Q 

E 

4- 

TO  P 

< 

OQ 

o 

OQ 

o 

o 

< 

< 

< 

1 

OQ 

AINi2Hf 

cF16 

(BiF3) 

0 

32 

24 

24 

AI3Ni12Hf 

cP4 

(AuCu3) 

0 

9 

0 

12 

AINi2Nb 

cF16 

(BiF3) 

0 

32 

24 

24 

AINi2Ti 

cF16 

(BiF3) 

0 

32 

24 

24 

AINiY 

hP9 

(Fe2P) 

3 

12 

18 

0 

AI4NiY 

oC24 

(AUNiY) 

52 

28 

48 

0 

AINi8Y3 

hP24 

(CeNi3) 

0 

18 

6 

48 

AI2Ni6Y3 

cl44 

(Ag8Ca3) 

0 

48 

48 

72 

b  Gas  standard  state. 
c  Metal  standard  state. 


0 

0 

■D 

0 

■D 

~o 

o 

o 

0 

0 

3 

O 

o 

'c 

'c 

CO 

TO 

ID 

3 

3 

0 

0 

0 

4_ 

0 

4- 

0 

E 

4_ 

a 

4_ 

a 

Q. 

o. 

CO 

~D 

CO 

~o 

N 

O  ~ 

>  ■ 

N 

>»  ■ 

N 

>  • 

C 

c 

GO  a) 

cq  a> 

CQ  0 

o 

n 

o 

1 

o 

-Q 

o 

J  o 
i  E 

I  -3 

o 

i  E 

I  -> 

J  o 
i  E 

I  -) 

OQ 

o 

_ <  i _ 

_ <  i _ 

_ <  £ 

32 

0 

-233±7 

-221 ±7 

-301 ±30 

3 

0 

-634±21 

-537±16 

-652±64 

32 

0 

-154±4 

-181  ±5 

-204±20 

32 

0 

-224 

-228±7 

-228±22 

15 

6 

-162±3 

-239±7 

-220±22 

12 

4 

-362±12 

-389±12 

-360±35 

78 

10 

-455±30 

-535±16 

-524±51 

96 

24 

-534±17 

-513±15 

-550±54 
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Table  S6  Second-neighbour  contribution  to  8M 


Element 

Stacking  Fault  Energy 
(mJ  m'2) a 

■o 

c 

o 

tq 

b.  > 

'w3  i 

Ag 

16 

0.0024 

0.0049 

Al 

166 

0.0238 

0.0416 

Au 

32 

0.0047 

0.0074 

Cd 

175 

0.0311 

0.1610 

Cu 

45 

0.0051 

0.0088 

Ir 

300 

0.0400 

0.0346 

Mg 

125 

0.0231 

0.0908 

Ni 

125 

0.0143 

0.0193 

Pd 

180 

0.0262 

0.0402 

Pt 

322 

0.0448 

0.0459 

Rh 

750 

0.0942 

0.0981 

115 

0.0263 

0.0253 

Zn 

140 

0.0198 

0.0878 

Approved  for  public  release;  distribution  unlimited. 


35 


D(AI)  =  141pm  Atom  Fraction  Hf  D(Hf)  =  158  pm 


0.00  0.20  0.40  0.60  0.80  1.00 

D(AI)  =  141pm  Atom  Fraction  Zr  D(Zr)  =  160  pm 


D(AI)=  141pm  Atom  Fraction  Ni  D(Ni)=126pm 


D(AI)  =  141pm  Atom  Fraction  Pt  D(Ni)  =  139  pm 


Figure  SI  |  Fractions  of  A-A,  A-B,  and  B-B  bonds  for  selected  binary  systems  as  a  function  of  atom  fraction  of  element  B  (fB).  Element  B,  where 
sbb  <  saa,  is  on  the  right  hand  side  of  each  binary  system.  Dashed  lines  indicate  the  bond  concentrations  expected  for  an  ideal  solution  of  equal¬ 
sized  atoms.  Constituent  atom  sizes  40  are  included  in  each  figure — note  the  universal  skewing  of  A-B  bonds  towards  the  smaller  constituent. 
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